Zero- Temperature Magnetic Transition in an Easy-Axis Kondo Lattice Model 
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We address the quantum transition of a spin- 1/2 antiferromagnetic Kondo lattice model with 
an easy-axis anisotropy using the extended dynamical mean field theory. We derive results in real 
frequency using the bosonic numerical renormalization group (bNRG) method and compare them 
with Quantum Monte Carlo results in Matsubara frequency. The bNRG results show a logarithmic 
divergence in the critical local spin susceptibility, signaling a destruction of Kondo screening. The 
T = transition is consistent with being second order. The bNRG results also display some subtle 
features; we identify their origin and suggest means for further microscopic studies. 
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A sizable number of (nearly) stoichiometric heavy 
fermions have recently been discovered in which the an- 
tiferromagnetic transition temperature can be continu- 
ously suppressed to zero [H. These materials have not 
only allowed further elucidation of the heavy fcrmion 
physics but also provided a concrete setting to address 
the larger question of quantum criticality. The applica- 
tion of the Landau paradigm considers the fluctuations 
of the magnetic order parameter as the primary criti- 
cal modes The resulting T = spin-density- wave 
(SDW) quantum critical point (QCP) @, % is Gaus- 
sian. However, a host of dynamical, transport and ther- 
modynamic data [H, H, @, 0, 0| suggest that the observed 
QCPs are non-Gaussian, indicating the existence of ad- 
ditional quantum critical modes. Since there is not yet a 
universal prescription for the identification of such emer- 
gent critical modes, microscopic considerations have been 
playing an important role. 

One idea invokes the breakdown of the Kondo screen- 
ing effect at the magnetic QGP to characterize the new 
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critical modes [13, lUj- In local quantum critical- 
ity Q , the destruction of the Kondo effect arises through 
the decoherence by the magnetic order parameter fluc- 
tuations. Microscopically, this picture has been stud- 
ied through the extended dynamical mean field theory 
(EDMFT) approach Here, the Kondo lattice sys- 

tems are analyzed in terms of a Bose- Fermi Kondo (BFK) 
model, with the spectra of its fermionic and bosonic 
baths self-consistently determined. The EDMFT ap- 
proach addresses the RKKY-Kondo competition, going 



beyond the seminal works of Refs. [IJ, ll5| in ways that 
are important for the collapse of the Kondo scale at the 
magnetic QCP. It treats this competition dynamically. 
Equally important, it incorporates not only paramag- 
netic/antiferromagnetic phases with a "large" Fermi sur- 
face, but also an antiferromagnetic phase with a "small" 
Fermi surface (local moments not participating in the 
electronic Fermi volume). The critical behavior of the 
BFK model was shown to allow [Si a self-consistent so- 



lution in which the criticality of the BFK model - with 
critical Kondo screening ~ is manifested at the magnetic 
QCP of the lattice. This analytical result was verified 
in a Quantum Monte Carlo (QMC) study of a Kondo 
lattice model with an easy-axis anisotropy [lij . An im- 
portant question is whether the actual zero-temperature 
transition is second order. Earlier works at finite temper- 
atures, using various QMC ap proaches, have led to some 
conflicting conclusions 17|,ll8|. The differences have been 
attributed to the different EDMFT equations, which han- 
dle the generated RKKY interactions on the ordered side 
differently 1^, 



In this Letter, we study the magnetic transition of the 
anisotropic Kondo lattice model directly at zero tempera- 
ture, using the recently developed bNRG method j2l].[2^. 
Our results are important for experiments, not only be- 
cause the numerical studies play an important role in 
the understanding of the unusual magnetic dynamics Q 
(which itself was the primary initial experimental indica- 
tion for the non-SDW nature of the QCP), but also be- 
cause the theoretical picture has crucial predictions for 
other experiments that are actively being examined by 
on-going experiments {e.g., Refs. [3, [1]). More gener- 
ally, whether unconventional QCPs would be stable and 
relevant to realistic models/materials or tend to be pre- 
emptied by first order transitions arc broadly important 
and also arises 12311 in, e.g., the case of dcconfined quan- 
tum criticality |24| in spin/boson lattice systems. 

The Kondo lattice Hamiltonian is 



Jk S,.s,,,+V(/,,/2) 5f5? (1) 



Here, St and Sc,i represent the spins of the S = ^ local 
moment and conduction c-electrons respectively. There 
are 1 moment and, on average, a; < 1 conduction elec- 
trons, per site. Jk is the antiferromagnetic Kondo in- 
teraction, tij is the hopping integral, corresponding to a 
band dispersion ek whose density of states (DOS) /3o(e) 
is featureless. In denotes the RKKY interaction; its 



Fourier transform, Jq, is the most negative at an antifer- 
romagnetic (AF) wavevector Q (/q = -/). The EDMFT 
approach leads to the effective impurity action [17 1 
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where iStop is the Berry phase of the local moment, and 
^loci G^^, and ^'^e the static and dynamical Weiss 
fields satisfying the self-consistency condition: 



h\oc = 

Xloc(w) = 
Gioc,(t(w) = 
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Here, M{lo) and Scr(w) are respectively the spin 
and conduction-electron self-energies, which satisfy the 
Dyson(-likc) equations: Sc,(w) = Go i.(a') - Gj^c.crl^)' 
and M{uj) = Xo^i'^) + Xioc('^)- '"af = (5^)™? is the 
staggered magnetization; xioc(w) and G\oc,rj{!^) are the 
connected local spin susceptibility and local conduction 
electron Green's function, respectively. Finally, M{lo) 
also specifies the lattice spin susceptibility [l^ : 



X(q,c^) = l/[/q + A/(^)] 



(4) 



As described in detail in Refs. , , the effective im- 
purity action [Eq. ^] can be rewritten in a Hamiltonian 
form, in which the dynamical Weiss fields arc represented 
by a fcrmionic bath and a bosonic one. Through a canoni- 
cal transformation, the fcrmionic coupling is reduced to a 
transverse field Ising model with an ohmic bosonic bath. 
Integrating out the two bosonic baths yields a form that 
is suitable for QMC studies: 



dT[hi,, 5--(T)+r5"(r) 
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Here, lCc(i^n) = ^da.'nl describes the ohmic dissipation; 
He and F are determined by the longitudinal and trans- 
verse components of the Kondo coupling, respectively. 

For the bNRG studies, we work in the real frequency 
domain by rewriting Eq. in a Hamiltonian form: 



'^imp ^loc S 
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(6) 




FIG. 1: (Color online) Inverse static AF susceptibility, Xaf = 
X(q = Q,a; = 0)"^ = M{uj = 0)-/, from the PM (trial /iioc = 
0; green circles) and AF (trial ftioc 7^ 0; blue circles) solutions 
and the AF order parameter ttiaf (red squares), obtained from 
the bNRG (a) and QMC (b) methods. The lines are guides 
to eyes. See the main text for details. 



where Cjp and gp are such that — S^-i^'^ = Xo ^('^) = 
Xo ^(w) - /Cc(w). The EDMFT procedure starts with a 
trial h\oc and Xo('^)- The bNRG iteration loop 2l| is 
then used to solve the impurity model © for rriAF and 
Xioc(^) which; in turn, lead to updated h\oc and Xo^i^)- 
The procedure is repeated until convergence is achieved. 
For the most part, we consider two-dimensional magnetic 
fluctuations [17| as represented by a constant RKKY 
DOS pi{e) = Eq 5{e - Q = (l/2/)e(/ - |e|), with 6 
being the Heaviside function. In this case, Eq. (|3b[) yields 



A/(cj) = //tanh[/xioc(w)] 



(7) 



We take the energy cutoff Wcutoff' = 1 and the parameters 
F = 0.75 and Kc = tt, yielding = l/xioc(t^ = 0,/ = 
0) « 0.71. In most cases (exceptions will be specified), 
we choose the NRG discretization parameter A = 2, keep 
Nt = 100 bosonic states for the impurity site and 8 states 
for the other sites, and retain = 60 many-body states. 
To reach convergence, the difference between two consec- 
utive iterations in hioc as well as in Xq ^(w) for each lu is 
smaller than 10~^. Away from the transition region, 30 
or so EDMFT iterations are sufficient. In the transition 
region, it takes as many as 2400 iterations. 

The resulting phase diagram is summarized in 
Fig. [Tfa). We observe a substantial drop of Xaf from 
both sides, as well as of ttiaf. Indeed, the magnetic order 
parameter ttiaf vanishes continuously within the numer- 
ical uncertainty as / approaches the transition point Id 
(« 1.1228r^). 

Fig. E^a) shows xioc(w) at various /, from around 
/ = Id and beyond. Above a cutoff scale, the local sus- 
ceptibility is logarithmically dependent on the frequency. 
Such a singular behavior signals the Kondo screening be- 
ing critical, which is the hallmark of local quantum criti- 
cality. Fitting the slope of the logarithmic dependence in 
terms of a/21 yields an a which is nearly constant (vary- 
ing by less than 2%) in the shown range of /. Through 
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FIG. 2: (Color online) Frequency dependence of the local 
spin susceptibility at various values of I around the magnetic 
transition, (a) Xiods^) vs. the real frequency u, from bNRG; 
(b) Xioc{ijJn) vs. the Matsubara frequency ujn, from QMC. 

the self-consistency Eq. ([7]) together with Eq. ([4]), a is 
the critical exponent that appears in the dynamical AF 
spin susceptibihty, Xaf(w) = x{Q,i^)- Extrapolating to 
the bNRG continuum limit (A 1+) yields a k, 0.83. 
The low-frequency cut-off scale for the logarithmic de- 
pendence is relatively small, becoming of the order of 
^ 10~^T]^ for the largest / we have reached as / is in- 
creased towards Ic2 (~ 1-26T]^), the instability point of 
the paramagnetic solution signaled by a diverging Xaf; 
this cut-off scale extrapolates to zero as / — *■ . 

For comparison, the QMC results for the phase dia- 
gram and the Matsubara frequency dependence of the lo- 
cal dynamical spin susceptibility arc shown in Figs, [ijb) 
and[2l^b), respectively. At the gross level, the bNRG and 
QMC results are similar to each other. 

At a fine level, the bNRG results contain some dif- 
ferences from their QMC counterparts. Major among 
these is the observation, as seen in Fig. [ija), that Ic2 
is larger than Id by about 12%. In the QMC re- 
sults, by contrast, Ic2 equals I^i within the numerical 
uncertainty of a few percents. To see whether this is 
unique to the bNRG results for the 2D magnetic fluc- 
tuations, we have carried out similar bNRG studies of 
the EDMFT phase diagram in the case of 3D magnetic 
fluctuations - as represented by a semicircular RKKY 
DOS, pi{e) = - e^e{I - |e|). The 3D case does 

not have the complication of a divergent local suscepti- 
bility, and a SDW solution is expected in the EDMFT 
approach \v^ . We find that the magnetic transi- 

tion is essentially continuous (with the upper bound of 
the order-parameter jump being 0.016), yet {Ic2^Ici) / Ici 
is still non-zero (about 13%). 

The observation of a continuous onset of the magnetic 
order parameter m^p but, at the same time, different 
Id and /c2, is unexpected. One possibility is that the 
dichotomy is inherent to the EDMFT equations. To ad- 
dress this, we return to the self-consistent equation for 
the magnetic order parameter, Eq. (|3ap . For a small h\oc 
(we have numerically determined that the magnetic so- 
lution is the same regardless of whether an infinitesimal 



finite value or a large value is chosen for the trial /iioc), 

we have = -Xioc^ioc - aa^Lc ~ (^^^loc + Note, 

the linear coefficient is equal to — Xioc(w = 0) since, in 
Eq. (O, /iioc couples linearly to only [g^]. We can 
then rewrite Eq. ([5a|) as 

rh\oc = -w/ifoc ~ ""^foc + ■■■ 1 (8) 

where r = XIoc/Xaf is the quadratic coefficient of the cor- 
responding static Landau function, and u = — asiXioc ~ 
Xaf] and v ~ ^ asiXioc ~ Xaf] are the quartic and sextic 
Landau coefficients. When m > (the alternative, u < 0, 
would lead to a large jump in ttiaf, in contrast to what 
we have observed) , we have a canonical case of a second- 
order transition at r = (in other words, a h\oc ^ solu- 
tion cannot occur for any r > 0). Through r ~ XIoc/Xaf, 
this implies that, at Id (the onset of the magnetic tran- 
sition) Xaf diverges. This is the same condition for 7^2, 
where the paramagnetic solution goes away. So, within 
the EDMFT equations per se, a continuous onset in TOaf 
must coincide with a vanishing {Ic2 — Ici ) ■ 

We are then led to search for numerical origins for the 
dichotomatic observation, and have identified the pri- 
mary source. Within bNRG, as in any NRG method, 
the imaginary part of the local susceptibility, x(oc('^)> 
calculated in terms of a set of Gaussian-broadened delta 
functions. The real part is in turn determined via the 
Kramers-Kronig relation, which we call x'loc kk('^)- T^^^ 
static local susceptibility can alternatively be calculated 
in terms of a) the differential response of the local mag- 
netization w.r.t. h\oc or b) X]„ l("l'5'z|0)P/(-B„ - £^0), 
where n labels all the many-body excited states and the 
ground state. We find that the latter two methods yield 
essentially the same result, which we call xioc, static- A key 
observation is that xioc.static is larger than Xioc,kk(^ = 0) 
by a sizable amount (about 11.5% for A=2, in the 2D 
case). The quadratic Landau coefficient then becomes 

r = Xloc,static/XAF - [Xloc,static/Xloc,KK(t^ = 0) - 1]. It fol- 

lows that the onset of the magnetic transition (at r ~ 0) 
already occurs before Xaf diverges, which explains the 
Id < Ic2 discussed earlier. In order to confirm our obser- 
vation, we have implemented the simplest modification 
scheme to ensure that the Kramers-Kronig of the NRG- 
calculated x(oc('^) yields a static local susceptibility that 
is equal to xioc, static- We use, during each EDMFT iter- 
ation, xioc, static for the a; = component of Xiocl^); 
retain x[oc kk('^) fo^ finite frequencies. We find that 
Id is increased compared to that of the "vanilla" scheme. 
(/c2 is essentially unchanged, although the normalization 
parameter T'^ is reduced.) Moreover, as shown in Fig. [3] 
for the 2D case, Ic2 ~ ^ci [with a difference less than 1% 
(2%) in the 2D (3D) case, with A = 2]. itiaf vs. {I -Id) 
from the modified scheme is mostly comparable to that 
of Fig. [Ija) except for being steeper in the immediate 
vicinity of Id [when (/ — Id) I Id is within a few per- 
cents] . The magnetic transition is therefore second order 
within the numerical accuracy. 
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FIG. 3: (Color online) Inverse AF susceptibilities (PM, green 
circles; AF, blue circles) and the AF ordered moment (red 
squares), from the modified scheme, as described in the main 
text. The notations are the same as in Fig. [T] 



For / > Ic2 l25|, the nominally self-consistent param- 
agnetic solution has max[2/x[oc('^)] > tt, which, through 
Eq. (H), yields an oscillatory M"{uj) [ig. By contrast, 
in the Matsubara frequency domain, a nominally para- 
magnetic solution (/iioc = = mAF, but a finite Curie 
constant) exists for / > Ic2, which helped to determine 
the phase diagram [l7| . 

Independently, Glossop and Ingersent [H] have car- 
ried out NRG studies within the EDMFT approach to 
the same Kondo lattice model. They used the NRG 
method of Ref. in which the Kondo coupling to 

the conduction electrons are directly treated (instead of 
being mapped to an Ohmic dissipation). Moreover, in 
Ref. (H^lb), they adopted a somewhat different modi- 
fication scheme to ensure the consistency between the 
static local susceptibilities from two ways of calculation 
within NRG. In spite of these differences in methods, the 
results from the two groups are largely compatible with 
each other. 

To summarize, we have carried out bosonic numerical 
renormalization group studies of the extended dynam- 
ical mean field theory of a Kondo lattice model. The 
local spin susceptibility has a logarithmic frequency de- 
pendence - signifying the critical Kondo screening - and 
the magnetic transition is consistent with being second 
order. These results provide evidence for local quantum 
criticality. Our study has also advanced the understand- 
ing of the numerical renormalization group, a venerable 
method [2^ in the area of correlated systems. 
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